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7 Abstract — A comparison of snow depths on sea ice was made 

8 using airborne altimeters and an Advanced Microwave Scanning 

9 Radiometer for the Earth Observing System (AMSR-E) simulator. 

10 The data were collected during the March 2006 National Aero- 

11 nautics and Space Administration (NASA) Arctic field campaign 

12 utilizing the NASA P-3B aircraft. The campaign consisted of an 

13 initial series of coordinated surface and aircraft measurements 

14 over Elson Lagoon, Alaska and adjacent seas followed by a se- 

15 ries of large-scale (100 km x 50 km) coordinated aircraft and 

16 AMSR-E snow depth measurements over portions of the Chukchi 

17 and Beaufort seas. This paper focuses on the latter part of the 

18 campaign. The P-3B aircraft carried the University of Colorado 

19 Polarimetric Scanning Radiometer (PSR-A), the NASA Wallops 

20 Airborne Topographic Mapper (ATM) lidar altimeter, and the 
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University of Kansas Delay-Doppler (D2P) radar altimeter. The 21 
PSR-A was used as an AMSR-E simulator, whereas the ATM and 22 
D2P altimeters were used in combination to provide an indepen- 23 
dent estimate of snow depth. Results of a comparison between the 24 
altimeter-derived snow depths and the equivalent AMSR-E snow 25 
depths using PSR-A brightness temperatures calibrated relative 26 
to AMSR-E are presented. Data collected over a frozen coastal 27 
polynya were used to intercalibrate the ATM and D2P altimeters 28 
before estimating an altimeter snow depth. Results show that the 29 
mean difference between the PSR and altimeter snow depths is 30 
—2.4 cm (PSR minus altimeter) with a standard deviation of 31 
7.7 cm. The RMS difference is 8.0 cm. The overall correlation 32 
between the two snow depth data sets is 0.59. 33 

Index Terms — Author, please supply index terms/keywords 34 
for your paper. To download the IEEE Taxonomy go to 35 
http://www.ieee.org/documents/2009Taxonomy_vl01.pdf. 36 AQl 

I. Introduction 37 

T he primary objective of the National Aeronautics and 38 
Space Administration (NASA) March 2006 Arctic field 39 
campaign was to assess the accuracy of the Aqua Advanced 40 
Microwave Scanning Radiometer for the Earth Observing Sys- 41 
tern (EOS) (AMSR-E) snow depth on sea ice retrievals [1]. 42 
The field campaign consisted of an initial series of coordinated 43 
surface and NASA P-3B aircraft measurements over Elson 44 
Lagoon, Alaska and adjacent seas on March 18 and 20 followed 45 
by a series of large-scale (100 km x 50 km) coordinated 46 
aircraft and Aqua AMSR-E measurements over portions of 47 
the Chukchi Sea, Kotzebue Sound, and the Beaufort Sea on 48 
March 21, 22, and 25, respectively. A sixth flight on March 24 49 
was coordinated with an ICESat overpass in the high Arctic 50 
to support a study of the effects of snow cover variability 51 
on ice thickness retrievals from the ICESat laser altimeter 52 
[2]. All six flights were made from Eairbanks International 53 
Airport, Alaska [Eig. 1(a)]. A transit flight to Greenland was 54 
also made on March 27 in coordination with an Envisat Radar 55 
Altimeter-2 overpass in the high Arctic to validate sea ice 56 
elevation measurements derived from the Envisat microwave 57 
altimeter [3]. 58 

The Elson Lagoon flights on March 18 and 20 were used 59 
to compare in-situ snow depth measurements with snow depth 60 
measurements made from the airborne radiometer and altime- 61 
ters. The results from these flights will be the subject of a 62 
forthcoming paper. In this paper, we use data collected over 63 
the flight areas of March 21, 22, and 25 [Eig. 1(a)] to compare 64 
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(b) 


Fig. I. (a) Six NASA P-3B flights made from Fairbanks, AK covered portions of Elson Lagoon near Pt. Barrow, AK, the Chucki and Beaufort seas, Kotzebue 4/C 
Sound, and the high Arctic during the March 2006 AMSR-E Arctic held campaign, (b) AMSR-E snow depth map (5-day average) for March 21, 2006. The color 
scale gives the snow depth in centimeters. Multiyear sea ice is masked out, because the snow depth retrievals are limited to first-year sea ice types only. 


65 the snow depth retrievals obtained from the NASA P-3B altime- 

66 ters and from the radiometer which has the same radiometric 

67 channels as the AMSR-E sensor. Even with the aircraft making 

68 two or three passes over an AMSR-E 12.5 km grid cell, the 


coverage by the aircraft sensors was too sparse for a direct com- 69 
parison with AMSR-E snow depths. Thus, we use the airborne 70 
radiometer as an AMSR-E simulator to compare the microwave 71 
radiometer and altimeter snow depths. Previous work used 72 
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TABLE I 

NASA P-3B Aircraft Sensors Flown During the Arctic 2006 Field Campaign 


Sensor 

Characteristics 

FOV at 680 ft 
Altitude 

Estimated 

Precision 

Purpose 

Polarimetric Scanning 
Radiometer (PSR-A 
and PSR-CX) 

Operating Frequencies 
(H&V-pol):6, 10, 18, 22, 
37, 89 GHz 

67 m 

1-2 K 

Aircraft AMSR-E simulator 
microwave based snow depth 
determination 

Airborne Topographic 
Mapper (ATM-II) 

Scanning Lidar altimeter 
combined with a 
differential GPS system 

2 m (100m 
swath width) 

10 cm 
[13] 

Air/snow interface elevation. 
Maps ice surface topography 
at high resolution 

Dual-Doppler Radar 
(Ku-band) 

Delay-Doppler Phase- 
monopulse (D2P) radar 
altimeter data 

30 m (across) x 
4 m (along) 

5 cm 
[12] 

Sea ice/snow interface 
elevation 

Digital Cameras 

2 KODAK Digital 
cameras (3 megapixel) 

170 m (across) 
X 260 m (along) 


Visible record of ice surface 


73 both the high-resolution airborne laser altimeter retrievals of 

74 snow-ice freeboard and the passive microwave retrievals of 

75 snow depth from this campaign to provide insight into the 

76 spatial variability of these quantities as well as optimal methods 

77 for combining high-resolution satellite altimeter measurements 

78 with low-resolution snow depth data [4]. 

79 The original intent of this work was to use the airborne 

80 altimeters as a validation tool to assess the AMSR-E sea snow 

81 on sea ice retrievals, but since the altimeter elevation differences 

82 used as a measure of snow depth on sea ice have yet to be vali- 

83 dated, we present a comparison between the airborne altimeter- 

84 derived snow depths and the airborne microwave radiometer- 

85 derived snow depths using an equivalent AMSR-E snow depth 

86 on sea ice algorithm. The comparative results provide insight 

87 into the limitations of both the altimetric and radiometric snow 

88 depth retrievals. 

89 II. Methodology 

90 A. EOS Aqua AMSR-E Satellite Data 

91 The AMSR-E was launched in May 2002 on the Aqua satel- 

92 lite. AMSR-E is a state-of-the-art sensor measuring microwave 

93 emissions over a broader range of wavelengths and with better 

94 spatial resolution than previous satellite radiometers. AMSR-E 

95 was designed and built by the Japan Aerospace Exploration 

96 Agency for the NASA EOS Aqua spacecraft [5]. The three 

97 AMSR-E sea ice products include sea ice concentration, snow 

98 depth on sea ice, and sea ice drift. In this paper, we make use of 

99 the snow depth on sea ice product. 

100 AMSR-E snow depth on sea ice is a 5 -day averaged gridded 

1 01 product at a resolution of 12.5 km and is derived using an 

102 algorithm described by [6]. While the product is available for 

103 both the Antarctic and Arctic, in the latter region, the snow 

104 depth retrievals are limited to areas of first-year sea ice, because 

105 multiyear ice presents a fundamental ambiguity, which is dis- 

106 cussed later, making the retrieval of snow depth over multiyear 

107 ice indeterminate, at least at present. An example of the 5 -day 

108 AMSR-E snow depth product is shown in Eig. 1(b). 

109 As described in [6], the snow depth on sea ice algorithm is 
no linearly related to the spectral gradient ratio corrected for sea 
III ice concentration GRV (ice) defined by 

GRV(ice) = [Tb(37V) - Tb(18V) - - C)] 

/ [Tb(37V) + Tb(18V) - k2{l - C)] (1) 


where Tb(37V) and Tb(18V) are the brightness temperatures 1 12 


of the satellite radiometer and 1 13 

ki = Tbow(37V) - Tbow(18V) (2) 

k2 =Tbow(37V) + Tbow(18V). (3) 

Tbow is the open water brightness temperature, and C is the sea 1 14 
ice concentration as determined by the enhanced NASA Team 1 15 
(NT2) algorithm applied to the AMSR-E data [7]. 1 16 

The snow depth in centimeters is given by 1 17 

/is = + U 2 GRV(ice). (4) 


Both the a\ and U 2 coefficients were derived from a lin- 1 18 
ear regression of in-situ snow depth measurements on SSM/I 1 19 
microwave measurements [6]. These coefficients were sub- 120 
sequently adjusted to take into account brightness tempera- I2I 
ture calibration differences between SSM/I and AMSR-E. Eor 122 
SSM/I equivalent GRV, a\ has the value of 2.9 cm, and U 2 has 123 
the value of —782 cm. 124 

The basis of the algorithm assumes that scattering increases 125 
with increasing snow depth and that the scattering efficiency is 126 
greater at 37 GHz than at 18 GHz. Eor snow-free first-year sea 127 
ice, the gradient ratio is close to zero, and it becomes more and 128 
more negative as the differential scattering increases resulting 129 
from an increase in snow depth and/or an increase in grain size. 130 
The upper limit for snow depth retrievals is about 50 cm which 13 1 
is a result of the limited penetration depth at 37 GHz [8]. 132 

The algorithm is applicable to dry snow conditions only. At 133 
the onset of melt, the emissivities of both the 18 GHz and the 134 
37 GHz channels approach unity (that of a blackbody) and 135 
the gradient ratio approaches zero initially before becoming 136 
positive. Thus, snow depth is indeterminate under wet snow 137 
conditions. Snow, which can be wet during the day, frequently 138 
refreezes during the night. This refreezing results in very large 139 
grain sizes, which results in a reduced emissivity at 37 GHz 140 
relative to 18 GHz, thereby decreasing GRV (ice) and thus I4I 
results in an overestimate of snow depth. These thaw-freeze 142 
events cause large temporal variations in the snow depth re- 143 
trievals. This temporal information is used in the algorithm to 144 
flag the snow depths as indeterminate from those periods with 145 
large fluctuations. As in-situ grain size measurements are even 146 
less frequently collected than snow depth measurements, the 147 
influence of grain size variations could not be incorporated 148 
into the algorithm. Because of diurnal melt-freeze cycles and 149 
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150 sporadic weather effects, AMSR-E daily snow depth products 

151 are 5 -day running averages. 

152 Because of the higher sensitivity of snow depth retrievals 

153 to ice concentrations less than 20%, the algorithm limits snow 

154 depth retrievals to ice concentrations between 20% and 100%. 

155 Ice concentrations less than 20% appear almost exclusively near 

156 the ice edge, so the total area excluded is relatively small. 

157 Both multiyear ice and deep snow on top of first-year ice 

158 result in increasingly negative values for the spectral GR [9]; 

159 therefore, the algorithm only retrieves snow depth in the sea- 

160 sonal sea ice zones. We currently use a dynamic GRV based 

161 filter which approximates the multiyear sea ice cover. This 

162 multiyear ice mask is defined on October 1 of each year as 

163 sea ice which has GRV values of less than —0.03. The same 

164 GRV test is done for each subsequent day, with the resulting 

165 classification being limited by the boundary of the previous 

166 day’s mask, with an allowance of a 1 pixel perimeter, to take 

167 into account the possible motion of the multiyear ice pack. 



168 B. Aircraft Data Sets 

169 The NASA P-3B aircraft carried the University of Col- 

170 orado Polarimetric Scanning Radiometer (PSR-A), the NASA 

171 Wallops Airborne Topographic Mapper (ATM) lidar altime- 

172 ter, and the University of Kansas Delay-Doppler (D2P) radar 

173 altimeter. The PSR-A was used as an AMSR-E simulator, 

174 whereas the ATM measured the range from the aircraft to the 

175 air/snow interface and the D2P measured the range from the air- 

176 craft to the sea ice/snow interface. The processing of the altime- 

177 ter measured ranges is quite complex and is discussed in detail 

178 elsewhere (e.g., [10]-[12]). The altimeter products used in this 

179 study are given as elevations measured in meters relative to a 

180 common geoid. The difference in altimeter elevations (ATM- 

181 D2P) was used to provide an independent estimate of snow 

182 depth. A summary of the aircraft instrument operating char- 

183 acteristics as well as the estimated precision of the altimeters 

184 obtained from previous field campaigns is presented in Table I. 

185 The method employed consisted of making three flights 

186 (March 21, 22, and 25) over large areas (100 km x 50 km) 

187 covering 32 AMSR-E grid elements (12.5 km on a side) on each 

188 day. The day before each of these flights, we utilized near real- 

189 time AMSR-E snow depth maps to plan the next day’s flight. 

190 On March 21, we covered an area in the Chukchi Sea which 

191 had a relatively shallow snow cover [Eig. 2(a)]. On March 22, 

192 we overflew an area in Kotzebue Sound which had a somewhat 

193 deeper snow cover [Eig. 2(b)], and on March 25, we flew over 

194 an area in the Beaufort Sea which had the largest apparent 

195 snow cover [Eig. 2(c)]. The orientation of each rectangular box 

196 in Eig. 2 matches the orientation of the flight lines shown in 

197 Eig. 1(a) for corresponding days. 

198 Eor the purpose of utilizing the PSR as an AMSR-E simula- 

199 tor, we calibrated the PSR 19 GHz V-pol. and 37 GHz V-pol. 

200 brightness temperatures relative to AMSR-E making use of all 

201 the data obtained for March 21, 22, and 25 resulting in a total 

202 of 96 data points (Eig. 3). The justification for using the PSR as 

203 a proxy for AMSR-E is the high correlation (0.94) between the 

204 AMSR-E and PSR GRV parameters (Pig. 4). 



(c) 


Fig. 2. AMSR-E snow depths for portions of (a) the Chukchi Sea overflown 
on March 21, (b) Kotzebue Sound overflown on March 22, and (c) the Beaufort 
Sea overflown on March 25. The red rectangle in each image indicates the ap- 
proximate area overflown by the NASA P-3B aircraft. Each rectangle measures 
4 by 8 12.5 km AMSR-E pixels. The color scale gives snow depths in cm. 

Once the PSR 19V and 37V brightness temperatures were 205 
converted to equivalent AMSR-E brightness temperatures 206 
using the regression equations shown in Pig. 3, the AMSR-E 207 
snow depth algorithm was applied [(1) and (4)] to obtain PSR 208 
snow depths. 209 

Pield airborne laser and radar altimeter measurements show 210 
that the difference between the ATM elevation and the D2P 211 
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Fig. 3. AMSR-E versus PSR regression plot for TB(I9V) (left) and TB(37V) 
(right). 




Fig. 5. NASA P-3 flight tracks (gray thin lines) on an Aqua MODIS image 
of Kotzebue Sound for March 22. The aircraft altimeter data coverage is also 
shown (black heavy lines). The segment highlighted within the large area 
of grey ice (segment A) off the Alaskan coast was used to determine the 
altimeter elevation statistics and the resulting offset between the ATM and D2P 
elevations. The inset images are captured from the onboard digital camera and 
show the character of the ice surface within the coastal polynya. Segment B is 
the portion of the flight track used for the profiles in Fig. 6. 


PSR GRV (37/19) 

Fig. 4. AMSR-E versus PSR GRV regression plot. 

212 elevation provides a snow depth estimate consistent with cli- 

213 matologies [14], because the ATM measures the elevation of 

214 the air/snow interface and the D2P measures the elevation 

215 of the snow/ice interface both relative to a common geoid. 

216 Before using the altimeters as an alternate means of providing 

217 estimated snow depths, we needed to calibrate them relative to 

218 each other over some sea ice surface with a known snow depth. 

219 Newly frozen leads or polynyas provide such a surface. The 

220 rationale is that the ATM and D2P elevations should match over 

221 newly formed ice because there is only a minimal snow cover, 

222 if any at all. An analysis of ATM and D2P elevations measured 

223 over frozen leads and polynyas on all three days showed that the 

224 area with a minimum ATM-D2P elevation variance (2.41 cm) 

225 occurred over the frozen coastal polynya on March 22. 

226 The mean difference was —9.93 cm indicating that we needed a 

227 10 cm offset in the D2P elevations to obtain agreement between 

228 the two altimeters. While we cannot be sure that there was no 

229 snow cover, without this offset there were 122 negative snow 

230 depths obtained with a maximum negative value of —12 cm, 

231 whereas with the offset there were only 17 negative values the 

232 largest being —2 cm. 

233 Fig. 5 shows an Aqua MODIS image with the NASA P-3 

234 flight tracks superimposed for March 22, 2006. Segment A of 

235 the flight track over the coastal polynya was used to intercal- 

236 ibrate the two altimeters. The three aerial photographs shown 

237 as insets in Fig. 5 confirm that this segment was comprised of 

238 newly formed sea ice. Fig. 6 shows the effect of the 10-cm offset 

239 as applied to the D2P elevations which brings the ATM and D2P 

240 elevations into better agreement over frozen leads in a portion 

241 of the March 22 flight (segment B on Fig. 5). 



Distance (km) 

Fig. 6. Portion of the March 22 flight (segment B on Fig. 5) shows that a 10- 
cm offset applied to the D2P elevations brings the ATM and D2P elevations 
into better agreement over frozen leads. 


Finally, for the purpose of obtaining a geolocated airborne 242 
sensor data set, the D2P altimeter data were chosen as the 243 
reference location. The ATM elevation and PSR brightness 244 
temperature data were averaged over a 35 m diameter circle 245 
around each given valid D2P point. The 35-m data sets were 246 
smoothed either to a 1-km length scale or to the 12.5-km 247 
AMSR-E grid scale for the comparison studies discussed below. 248 


III. Results and Discussion 


249 


The sea ice and snow cover characteristics of the areas 250 
overflown on March 21, 22, and 25 are all quite different and are 251 
discussed in the context of their microwave polarization (PR) 252 
and spectral gradient (GR) signatures. PR is defined in terms of 253 
the 19-GHz horizontal and vertical polarization PSR channels 254 

PR19 = [TB19V - TB19H]/[TB19V + TB19H]. (5) 
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Fig. 7. Plots illustrate the differences in PSR microwave PR-GR signatures for 
the three study areas on (a) March 21, (b) March 22, and (c) March 25, 2006. 
In each plot, the locations of pure first-year (FY), new (NEW), and multiyear 
(MY) ice types are indicated. 

255 Whereas GR is defined in terms of the 19-GHz and 37-GHz 

256 vertical polarization PSR channels 

GRV37/19 = [TB37V - TB19V/[TB37V + TB19V]. (6) 

257 The PR-GR characteristics of each of these three areas are 

258 shown in Fig. 7 through the use of PR-GR scatter plots. The 

259 PR-GR plot for March 21 [Fig. 7(a)] shows a fairly tight cluster 

260 near PR of 0.05 and GRV of —0.02 which is typical of first-year 

261 ice types (e.g., [9]; [15]). A looser cluster of points, typical of 

262 new and young ice types, straddles the GRV value of 0 and 

263 extends to higher PR values. The plot for March 22 [Fig. 7(b)] 



■ 10-1 ^ ^ 

-0.1 -0.0S 0 0.05 

PSR GRV37/19 


0.1 


(b) March 22 Kotzebue Sound flight 



-10-3 

-0,3 


-0,05 


PSRGRV37/19 
(c) March 25 Beaufort Sea flight 


0,05 


0.1 



- 5 -^ ^ ^ ^ 1 

-0.1 -0.05 0 0.05 0.1 

PSR GRV37/19 

Fig. 8. Plots illustrate the relationship between the altimeter measured snow 
depths and the PSR GRV signatures for the three study areas on (a) March 21, 

(b) March 22, and (c) March 25, 2006. 

shows that in addition to the typical first-year ice distribution 264 
of points, many points have more negative GRV values. The 265 
more negative GRV values are likely the result of deeper snow 266 
and the effects of the melt/freeze event that occurred in mid 267 
February which may have resulted in a snow cover with ice 268 
layers resulting in more scattering of the 37-GHz radiation 269 
relative to 19 GHz. Finally, the area overflown on March 25 270 
was comprised of first-year and multiyear sea ice with no new 271 
and young ice types [Fig. 7(c)]. 272 

Scatter plots of the altimeter snow depths versus the PSR 273 
GRV values for each of the three study areas overflown are 274 
shown in Fig. 8. The expected linear relationship between the 275 


CAVALIERI et al.\ COMPARISON OF SNOW DEPTH ON SEA ICE RETRIEVALS 


7 



Fig. 9. Sequence of images showing daily-averaged ECMWF surface atmospheric temperatures (top) and AMSR-E snow depth retrievals (bottom) for a two- 4/C 
week period in Eebruary 2006. The study areas overflown on March 21 and 22 are indicated by red rectangles. 


276 microwave parameter GRV, which is the independent variable 

277 in the snow depth algorithm [6], and the altimeter snow depth is 

278 lost for the March 22 and March 25 areas. Only for the March 

279 21 area does the linear relationship hold [Fig. 8(a)]. 

280 Reasons for the lack of correlation for March 22 and 25 

281 [Fig. 8(b) and (c)] are difficult to determine with certainty. 

282 The lack of correlation for the March 25 flight in the Beaufort 


Sea is probably related to the large fraction of multiyear ice in 283 
the region. However, the March 22 area in Kotzebue Sound is 284 
devoid of multiyear ice, but contains ice having more negative 285 
GRV values [Fig. 7(b)] than is normally observed in first-year 286 
ice regions. As noted earlier, there was a large-scale melt-freeze 287 
event in Kotzebue Sound during mid-February 2006. Fig. 9 288 
shows a sequence of daily ECMWF (ERA-interim) atmospheric 289 
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^5 - 0.045 

Jan Feb Mar 

2006 


Fig. 10. Time series of the first 90 days of 2006 showing 6-hourIy ECMWF surface air temperatures and daily AMSR-E GRV values corresponding to the 
highlighted pixels in Fig. 9 for March 21 and 22. February 15, 2006 is the day when the air temperature exceeded 0 C. Note the difference in the behavior of the 
AMSR-E GRV values for the two flight regions after the onset of melt. 

290 temperatures and AMSR-E snow depth maps during mid- 

291 February covering both the Chukchi Sea and Kotzebue Sound 

292 flight areas overflown on March 21 and 22, respectively. It is 

293 clear from Fig. 9 that the flight area over Kotzebue Sound 

294 had positive daily- averaged air temperatures in mid-February, 

295 whereas the flight area over the Chukchi Sea had not. In partic- 

296 ular, during the period from February 14-20, two low pressure 

297 systems initially centered over the Gulf of Anadyr [Fig. 1(a)] 

298 migrated into the Chukchi and Beaufort seas, resulting in a 

299 combination of southerly winds, increased air temperatures, 

300 and a likely increase in down-welling, long-wave radiation as- 

301 sociated with increased cloud cover. With air temperatures near 

302 zero, it is also possible that some precipitation may have fallen 

303 as rain, which would have signiflcantly affected the scattering 

304 properties of the snow cover. The Kotzebue weather station 

305 reported warming daily air temperatures from the beginning of 

306 February to February 15 when the maximum temperature of 

307 1.1 °C was reached. An increase in snow on the ground was not 

308 reported until February 25 when the measured snow and ice 

309 on the ground doubled to 28 cm. The maximum snow cover of 

310 43 cm reported at Kotzebue was reached during mid-March. 

311 The weather conditions and melt event in Kotzebue Sound 

312 may have resulted in a combination of deep snow and a 

313 metamorphosed snow cover with ice layers and coarser-grained 

314 snow. This melt event which affected the entire flight area is a 

315 probable cause for the lack of correlation shown in Fig. 8(b). 

316 Fig. 10 provides a time series of 6-hourly ECMWF surface 

317 air temperatures [16] and daily AMSR-E GRV values for the 

318 highlighted (red) pixels shown in Fig. 9 for the first three 

319 months of 2006. The red pixel within the red rectangle for 

320 Kotzebue Sound is located in the upper left portion of the 

321 flight area, and the red pixel for the Chukchi Sea is in the 

322 upper portion (Fig. 9). Following February 15, 2006, the day 

323 of maximum air temperature (+1.08 C), there is a marked 

324 difference in the behavior of the AMSR-E GRV values for the 

325 two flight regions after the onset of melt. The Chukchi Sea 


region apparently did not undergo the same degree of surface 326 
melt on February 15 (Fig. 10). In fact, none of the 32 grid cells 327 
overflown on March 21 had daily average air temperatures 328 
above —0.9 C with the warmest temperatures occurring closest 329 
to Kotzebue Sound [upper left in Fig. 2(a)]. The average of 330 
the daily air temperatures on February 15 for the 32 grid cells 331 
overflown on March 21 was — 1.4 C. The GRV values for both 332 
regions decreased initially after the melt event. The Chukchi 333 
Sea GRV values became less negative beginning on March 334 
12 and maintained values between —0.005 and —0.01 from 335 
March 14 through March 29 (Fig. 10). The GRV values in this 336 
range are typical of new, young, and thin first-year ice types. 337 
Because the Chukchi Sea region is much more dynamic than 338 
Kotzebue Sound, one possibility is that the February Chukchi 339 
Sea ice cover was displaced by sea ice having different 340 
(younger) surface characteristics. To explore this possibility, 341 
we compare daily AMSR-E snow depth maps with IFREMER 342 
(Institut Fran 9 ais de Recherche pour 1’ exploitation de la Mer, 343 
Issy-les-Moulineaux, France) AMSR-E sea ice drift maps ob- 344 
tained from (ftp://ftp.ifremer.fr/ifremer/cersat/products/gridded/ 345 
psi-drift/documentation/amsr.pdf) for a 10-day period in March 346 
2006. These maps are shown in Fig. 11. 347 

From March 13 through March 17 the sea ice drift was 348 
toward the north, but from March 18, 19, and 20, there was 349 
even stronger ice drift away from the Alaskan coast (Fig. 11). 350 
The Alaskan coast region between Cape Lisbume and Point 351 
Lay [Fig. 1(a)] produces a large volume of ice each winter 352 
through oceanic heat loss by coastal polynyas. The ice produced 353 
is often swept up in large-scale cyclonic or anticyclonic gyres 354 
and transported to other parts of the Arctic Ocean. The snow 355 
depth maps in Fig. 1 1 show an increasingly large area of ice 356 
with a shallow snow cover. Presumably, recently formed new 357 
and young ice types were advected into the area overflown on 358 
March 21 resulting in less negative GRV values (Fig. 10). 359 

Next, we examine the AMSR-E pixel-averaged D2P and 360 
ATM elevations, the altimeter and PSR snow depths, the 361 
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Fig. 1 1 . Sequence of images showing IFREMER AMSR-E sea ice drifts for a 2-day period together with the AMSR-E snow depths from March 13 to March 22 4/C 

in the vicinity of Kotzebue Sound and the Chukchi Sea. The overflight areas for the Chukchi Sea on March 21 and for Kotzebue Sound on March 22 are indicated 
by red rectangles as in Fig. 9. 


362 ATM-derived surface roughness, and the AMSR-E snow depths 

363 for both the Chukchi Sea region overflown on March 21 

364 (Table II) and the Kotzebue Sound region overflown on March 

365 22 (Table III). The orientation of the AMSR-E grid elements 

366 in Table II is rotated 90° relative to the AMSR-E cells shown 

367 in Eig. 2(a). The orientation of the grid elements in Table III is 

368 similar to that shown in Eig. 2(b). The surface roughness was 


obtained by calculating the average standard deviation of the 369 
ATM elevations over each AMSR-E grid cell in each table. 370 
In Table II, for the Chukchi Sea area, both the D2P and ATM 371 
elevations show similar spatial patterns as do the altimeter and 372 
PSR snow depths with the deepest snow found in the upper left 373 
and lower right portions of the 32-cell grid. A comparison of 374 
the ATM roughness values with the altimeter and PSR snow 375 
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TABLE II 

Mean (A) D2P Elevation, (B) ATM Elevation, (C) Altimeter Snow Depth, (D) PSR Snow Depth, (E) ATM Roughness, and (F) AMSR-E 
Snow Depth eor Each op the 32 AMSR-E Grid Elements (Column, Row) Overplown on March 21, 2006. Shades op Gray From Light 
TO Dark are Used to Indicate Increasing Values From Low to High. There was no Aircrapt Coverage op Grid (377,156) 


AJ D2P Elevation 

370 

371 

372 

373 

374 

375 

376 

377 

154 

0.6091 

0.6550 

0.4342 

0.4272 

0 4422 

0 3525 

0.3446 

0.2287 

155 

0.4924 

0.4986 

0.5455 

0.6505 

0,4602 

0.4817 

0,4313 

0.4196 

156 

0.4353 

0.6096 

0.7945 

0.7141 

0 6537 

0 5595 

0.5269 


157 

0,4060 

0.6141 

0.7147 

0.6279 

0.6293 

0,66891 

0,5432 

0,5194 

B)ATM Elevation 

370 

371 

372 

373 

374 

375 

376 

377 

154 

0.9657 

0.8950 

0.5507 

0.5879 

0 6530 

0 5395 

0.5522 

0.4453 

155 

0.6550 

0,6419 

0.7050 

0,8923 

0.6657 

0.6488 

0.6287 

0.6806 

156 

0.6669 

0.8326 

1.0683 

0.9583 

0.8632 

0.7195 

0,7236 


157 

0.5284 

0.8244 

0.9221 

0.7715 

0 7593 

0.8433 

0,7359 

0.8126! 

C) Altimeter SD 

370 

371 

372 

373 

374 

375 

376 

377 

154 

35.66 

24.00 

11.65 

16.07 

21.09 

18,70 

20.76 

21.66 

155 

16.25 

14.33 

15.95 

2417 

18.54 

16.72 

19.73 

26.09 

156 

23.15 

22.30 

27.39 

24,43 

20.95 

16.00 

19.67 


157 

12.24 

21.03 

20,74 

14.36 

13.00 

17.44 

19.27 

29.32 

D) PSR SD 

370 

371 

372 

373 

374 

375 

376 

377 

154 

24.83 

18.97 

9.91 

12.78 

16.78 

14.25 

17 18 

19.33 

155 

11,27 

9.41 

11.91 

16 39 

13.62 

14.21 

18.57 

25.17 

156 

16.27 

18.73 

21,33 

19.87 

20.34 

18,43 

19.22 


157 

11.99 

20.70 

22.67 

16 98 

15.19 

17.94 

19.01 

26.27 

E)ATM Roughness 

370 

371 

372 

373 

374 

375 

376 

377 

154 

21.86 

22.21 

9.71 

16.68 

16.60 

18.18 

19.53 

17.27 

155 

17.08 

14.70 

16.38 

25.15 

17.81 

18,08 

18.19 

23.89 

156 

19.21 

22.71 

25.09 

24 27 

21.07 

15.58 

18.00 


157 

16.21 

26.64 

25,86 

20.85 

17.78 

16,85 

20.37 

21.65 1 

F) AMSR-E SD 

370 

371 

372 

373 

374 

375 

376 

377 

154 

19.40 

16.86 

13 42 

14 28 

16,24 

16.91 

19.00 

19.42 

155 

16.40 

13.76 

12,73 

15.07 

17.56 

15.73 

20.01 

25.71 

156 

14.72 

17.57 

19.27 

18 39 

19.20 

18.51 

17.08 


157 

16,23 

22.41 

25.03 

19.30 

18.73 

20.36 1 

19.40 

21,35 1 


TABLE III 

Mean (A) D2P Elevation, (B) ATM Elevation, (C) Altimeter Snow Depth, (D) PSR Snow Depth, (E) ATM Roughness, 
AND (F) AMSR-E Snow Depth eor Each op the 32 AMSR-E Grid Elements (Column, Row) Overplown on March 22, 2006. 
Shades op Gray From Light to Dark are Used to Indicate Increasing Values From Low to High 


A) D2P Elevation 

136 

137 

138 

139 

140 

141 

142 

143 

364 

1.810 

1.664 

1.436 

1.071 

0.704 

0.454 

0.306 

0.300 

365 

1.773 

1.552 

1.209 

0,850 

0.566 

0.361 

0.281 

0.240 

366 

1.719 

1.419 

1 012 

0.661 

0.510 

0 397 

0.251 

0.300 

367 

1.719 

1.291 

0.931 

0,646 

0.542 

0.468 

0.323 

0.361 

B)ATM Elevation 

136 

137 

138 

139 

140 

141 

142 

143 

364 

1.928 

1.752 

1 546 

1.205 

0.869 

0 594 

0.453 

0.485 

365 

1.880 

1.665 

1.328 

0.951 

0.701 

0.525 

0.438 

0.412 

366 

1.857 

1.559 

1.140 

0.806 

0.664 

0 570 

0.382 

0.484 

367 

1.893 

1.388 

1.087 

0.833 

0.716 

0.638 

0.470 

0.457 

CJAItimeter SD 

136 

137 

138 

139 

140 

141 

142 

143 

364 

11.855 

8.813 

10.973 

13.407 

16.456 

13.988 

14.683 

18,442 

365 

10.686 

11.256 

11.883 

10.115 

13.502 

16.417 

15.712 

17,179 

366 

13.805 

14.020 

12.791 

14.524 

15 373 

17.243 

13.119 

18.391 

367 

17.415 

9.721 

15.519 

18.641 

17.456i 

17.057 

14.688 

9.612 

D) PSR SD 

136 

137 

138 

139 

140 

141 

142 

143 

364 

43.2 

47.6 

29.8 

36,0 

30.6 : 

20,7 

19.6 

24.6 

365 

39,7 

40.4 

34 7 

34.9 

26.3 

18 7 

24.8 

28 3 

366 

27.4 

27.5 

24.3 

23.0 

22.2 

24.1 

22.3 

23.0 

367 

20.2 

20.2 

21.5 

23,7 

23.3 

28,5 

22.5 

14.8 

E)ATM Roughness 

136 

137 

138 

139 

140 

141 

142 

143 

364 

15.92 

13,29 

17.69 

21.19 

25,59 

22.48 

23,79 

26.51 

365 

15.02 

14 30 

15.62 

16 63 

22.48 

23.91 

24.50 

23.49 

366 

19.41 

19.45 

19.56 

22.10 

24.14 

24.90 

20.54 

26.62 

367 

25.60 

17,73 

23.13 

27.23 

26,19 

23.90 

22,37 

19.21 

F) AMSR-E SD 

136 

137 

138 

139 

140 

141 

142 

143 

364 

36.44 

39.02 

30.32 

33.32 

28,04 

22.69 

22,73 

24.27 

365 

35.77 

36.00 

32.60 

30.55 

26.24 

22.66 

24.57 

25.83 

366 

20.08 

26,00 

24.67 

24.57 

23,19 

22.96 

23,97 

20.45 

367 

19.25 

22,56 

21,61 

20.87 

21,08 

22.33 

18,57 

15.11 
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376 depths shows that there is a positive correlation between snow 

377 depth and surface roughness for both the altimeter and PSR 

378 distributions. This is consistent with previous studies (e.g., 

379 [17]). The AMSR-E snow depths are only weakly correlated 

380 with the surface roughness and the altimeter and psr snow 

381 depths. The latter result is probably due to the spatial sampling 

382 difference between aircraft and spacecraft. 

383 In Table III, for the Kotzebue Sound area, both the D2P 

384 and ATM elevations show a similar pattern with an increase 

385 in elevation from right to left which probably corresponds to a 

386 changing geoid. The change is about 1.5 m over a distance of 

387 1 00 km, length of the P-3 flight line (eight 12. 5 -km AMSR-E 

388 grid cells). A comparison of the altimeter and PSR snow depths 

389 shows no agreement for this particular day. In fact, there is 

390 deeper snow derived from the altimeters on the right side of the 

391 flight area, whereas the PSR deep snow is found on the left side 

392 of the area. One possible explanation is that the greatest effects 

393 from the mid-Eebruary melt/freeze event and storm passages 

394 were felt in the upper left of the flight area (see Eig. 9). Because 

395 of this large-scale event, the sea ice snow cover in the upper 

396 left portion of the flight area may have had ice layers imbedded 

397 in the snow cover, which would have been particularly likely 

398 if rainfall had occurred. These ice layers may have resulted 

399 in lower altimeter snow depths (Table III). Larger size snow 

400 grains in the affected area would have also caused the PSR snow 

401 depths to be overestimated [18], because of greater scattering at 

402 37 GHz relative to 19 GHz. Unfortunately, we do not have in- 

403 situ measurements to confirm this interpretation. 

404 Another factor influencing the altimeter snow depth retrievals 

405 is the change in velocity of electromagnetic radiation from air 

406 to snow. The snow depth correction (v/c), where v is the wave 

407 velocity in snow, c the speed of light in vacuo, is proportional to 

408 a/ {s'), where s' is the dielectric permittivity of saline snow (i.e., 

409 the real part of the dielectric constant). A dielectric mixture 

410 model for saline snow [19] has been used to compute s' . The 

41 1 model parameterization is a function of snow properties (den- 

412 sity p, salinity S, and temperature T), and the frequency of the 

413 radiation (15 GHz in our case). Our v/c correction ranges be- 

414 tween 0.7 (p = 400 kg/m3 5' = 15 ppt T = 265 K) and 0.8 (p = 

415 300 kg/m3 *5 = 0 ppt T = 255 K). This range has been used to 

416 establish uncertainties of the altimeter snow depths (Pig. 12). 

417 We plot the PSR snow depths versus the altimeter snow 

418 depths in Pig. 12 for the Chukchi Sea flight on March 21 

419 where we have a total of 880 coincident altimeter and PSR 

420 measurements spanning portions of 31 AMSR-E pixels. Por 

421 the purpose of gaining insight into the effects of the air/snow 

422 velocity differences on the snow depth retrievals, we show three 

423 regression lines, one for the uncorrected altimeter snow depths 

424 (dashed line) and two others for the corrected altimeter snow 

425 depths (using the 0.8 and 0.7 v/c factors). The uncorrected 

426 velocity has the smallest slope of 0.43, whereas the 0.7 and 0.8 

427 corrected retrievals have slopes of 0.54 and 0.62, respectively. 

428 Although these corrections increase the slope slightly, we still 

429 have slopes much less than 1 . The length of the error bar for 

430 each point shown in Pig. 12 is determined from the 0.7 and 

431 0.8 v/c corrections and provides a sense of how much the 

432 correction affects the snow depth retrieval. The variation in 

433 v/c which depends on the snow properties certainly contributes 


to the observed scatter. We also indicate surface roughness, 434 
which is computed from ATM measurements, for each data 435 
point in Pig. 12 through the use of a color scale. It is apparent 436 
that both the PSR and altimeter snow depths increase with 437 
increasing surface roughness. The correlations between the 438 
PSR and altimeter snow depths and surface roughness are 0.60 439 
and 0.67, respectively. 440 

Pinally, we calculate comparison statistics based on the PSR 441 
and altimeter snow depth data sets for the Chukchi Sea flight on 442 
March 21. We have not corrected the altimeter snow depths for 443 
air/snow velocity changes, because of the large uncertainty in 444 
the snow parameters needed for the correction. These statistics 445 
are presented in Table IV. The mean snow depth difference 446 
(PSR minus altimeter) is —2.4 cm with a standard deviation 447 
of 7.7 cm. The RMS error is 8.0 cm, and the overall correlation 448 
between the two snow depth data sets is 0.59. 449 

IV. Summary and Conclusions 450 

Although the original intent of the Arctic 2006 field cam- 451 
paign was to use the airborne altimeters as a validation tool to 452 
assess the AMSR-E snow on sea ice retrievals, we could not 453 
undertake a validation study, because the altimeter elevation 454 
differences as a measure of snow depth on sea ice have yet to 455 
be validated. Thus, we could not justifiably use the altimeter 456 
snow depths as a validation data set. Nonetheless, a com- 457 
parison between the altimeter-derived and radiometer-derived 458 
snow depths provided insight into the limitations of both 459 
approaches. 460 

Of the three flights made over the ice-covered seas surround- 461 
ing Alaska, only the Chukchi Sea flight made on March 21 462 
provided data which yielded a good correlation between the 463 
altimeter and radiometer snow depths. However, the slope of 464 
the regression line is much less (~0.5) than 1. An understanding 465 
of this requires a careful comparison of both the altimetric 466 
and radiometric retrieval methods with in-situ snow depth 467 
measurements. Snow depth retrievals over Kotzebue Sound on 468 
March 22 were apparently affected by a melt-freeze event in the 469 
previous month. This event may have produced ice layers in the 470 
snow cover resulting in an underestimate of snow depth by the 471 
altimeters. The first two flights were over first-year ice, whereas 472 
the third flight over the Beaufort Sea on March 25 covered 473 
an area comprised mostly of multiyear ice. The presence of 474 
multiyear ice results in an ambiguous radiometric snow depth 475 
signature, because of scattering of the upwelling radiation by 476 
empty brine pockets in the freeboard layer of the multiyear ice 477 
[20] . It is this ambiguous signature that probably led to the poor 478 
correlation between the two snow depth data sets. Currently, 479 
there is no way to distinguish between first-year ice with a deep 480 
snow cover and multiyear ice. 481 

The potential to retrieve snow depth from airborne lidar 482 
and radar altimeter measurements has been demonstrated in 483 
several studies (e.g., [12], [14]), but a true validation of this 484 
method has not yet been demonstrated. Purthermore, there is 485 
a recurrent need to apply an adjustment to the radar altimeter 486 
data. Indeed, over some areas, the surface (i.e., the air/snow 487 
interface) elevation tracked by the lidar is lower than the 488 
snow/ice interface that should be detected by the radar, resulting 489 
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4/C Fig. 12. PSR snow depths versus the airborne altimeter-derived snow depths for March 21, 2006. There are three regression lines: one for the uncorrected 

altimeter snow depths (dashed line), one each for the 0.7 v/c corrected (light solid line), and the 0.8 v/c corrected (dark solid line) altimeter snow depths. ATM- 
derived surface roughness for each point is color coded. 


TABLE IV 

Comparison Snow Depth Statistics for the March 21, 2006 Chukchi Sea Study Area 


Snow Depth 
Data Set 

Snow Depth 

Mean Snow Depth 

Mean Diff. 

RMS 

Range 

+ 1 Stan. Dev. 

(PSR-Alt SD) 

Diff. 

(cm) 

(cm) 

+ 1 Stan. Dev. (cm) 

(cm) 

PSR (Equiv. AMSR-E) 

11.9-44.4 

17.3 + 6.9 

-2.4 + 7.7 

8.0 

Altimeter (ATM-D2P) 

1.9-53.7 

19.7 ±9.3 




490 in unrealistic negative snow depths. An explanation of these 

491 negative snow depths is problematic. The current study and 

492 several previous ones [12], [14], [21] have encountered the need 

493 to adjust the radar altimeter measurements relative to the lidar 


measurements. Understanding this recurrent discrepancy must 494 
be a priority for future studies that aim at using the difference 495 
between airborne lidar and radar altimeter measurements as a 496 
proxy for snow depth. 497 
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498 Finally, the current status of the AMSR-E snow depth algo- 

499 rithm validation is that it is incomplete. We cannot provide an 

500 overall estimate of accuracy with any confidence. Some valida- 

501 tion studies have been undertaken with in-situ and ship-borne 

502 measurements [8], [17], but comparisons between satellite re- 

503 trievals and surface point measurements can in itself introduce 

504 biases [6]. Thus, there is still a critical need to develop validated 

505 methods of retrieving snow depth from airborne sensors to 

506 help bridge the spatial divide between satellite observations 

507 and surface point measurements. Furthermore, the AMSR-E 

508 snow depth algorithm currently does not take into account 

509 surface roughness or snow grain size variations, even though 

510 both of these parameters affect snow depth retrievals [8], [17]. 

5 1 1 More comparative studies are needed covering different surface 

512 conditions at different times of the year. Previous studies [8], 

513 [18], [22] suggest that the use of the 10-GHz AMSR-E channels 

514 may help both in differentiating between smooth and rough 

515 surfaces and in lessening the affect of increasing snow grain 

516 size. Thus, work remains to be done to improve snow depth 

517 on sea ice retrievals from both microwave radiometers and 

518 altimeters. 
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